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Abstract 

A continuum extrapolation of static two- and four-quark energies calculated in 
quenched SU(2) lattice Monte Carlo is carried out based on Sommer's method 
of setting the scale. The /3-function is obtained as a side product of the ex- 
trapolations. Four-quark binding energies are found to be essentially constant at 
(3 > 2.35 unlike the two-body potentials. A model for four-quark energies, with 
explicit gluonic degrees of freedom removed, is fitted to these energies and the 
behaviour of the parameters of the model is investigated. An extension of the 
model using the first excited states of the two-body gluon field as additional basis 
states is found to be necessary for quarks at the corners of regular tetrahedra. 
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1 Introduction 



Systems of many hadrons play a crucial role in nature and it is important to 
understand the hadronic interactions from first principles. This would shed light 
on multiquark bound states and e.g. meson-meson scattering. A convenient tool 
for analysing the problem would be a potential model for multi-quark systems 
with explicit gluonic degrees of freedom removed. Such a model might be based 
on an effective two-body interaction as in the case of valence electrons in metals 
and nucleons in nuclei. 

Perturbation theory of QCD cannot even treat the confinement of quarks and 
gluons into hadrons, and the only known way to make realistic calculations of in- 
teracting quark clusters are Monte Carlo lattice techniques. Green and coworkers 
have been simulating systems of four static quarks in quenched SU(2) mainly on 
16 3 x 32 lattices references therein]. Four quarks, because it is the smallest 
number that can be partitioned into different colour-singlet groups. Energies of 
several configurations such as rectangular (R), linear (L) and tetrahedral (T) - 
examples of which are shown in figure |l] - have been simulated to get a set of 
geometries representative of the general case. The j3 values used have been 2.4 
and 2.5. 

However, e.g. Booth et al. conclude that asymptotic scaling for SU(2) gauge 
theory begins at (3 > 2.85, so there is reason to suspect that the simulation results 
could contain significant lattice artefacts. Previous estimates show that finite 
size effects of our simulations are unimportant, while the conclusion that scaling 
(as opposed to asymptotic scaling) has been achieved at (5 — 2.4 is somewhat 
more questionable for larger configurations and excited states ||. It is therefore 
expected that a continuum extrapolation could possibly yield different energies 
and point out artefacts in a parameterization of the energies while making the 
physical content clearer. Such an extrapolation is the object of this work. 

The geometries to be simulated here were chosen to be squares and tilted rectan- 
gles [(R) with x = y and (TR) in figure |J, because the simple model described 
below works for them and these geometries also exhibit the largest binding en- 
ergies. The P values (and lattice sizes) used were f3 = 2.35, 2.4 (16 3 x 32), 2.45 
(20 3 x 32), 2.5 (24 3 x 32) and 2.55 (26 3 x 32). In addition to these two-body runs 
were performed at (3 — 2.3 (12 3 x 32) to set the scale at a lower value of f3. 

2 A model for four-quark energies 

A model for the energy of four static quarks with explicit gluonic degrees of free- 
dom removed has been developed by Green and coworkers fl], references therein]. 
In this model a potential matrix is diagonalized with different two-body pairings 
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Figure 1: Some of the simulated four-quark geometries and their two-body pair- 
ings. 



as basis states. The basis states for some simulated configurations are shown in 
figure |l| For example, in the case of two basis states A and B, the eigenvalues A, 
are obtained from 

[V - AjN] % = 0, (1) 

with 

N=f \ IAT { /Nc ) and V=( V lt V2i ^ AB ), (2) 

where Vij represents the static two-body potential between quarks i and j. Vab 
comes from the perturbative expression 

= —Af(N c )Ti ■ T j% , (3) 
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where for a colour singlet state the normalization is chosen to give 

< >= v ij- The four-quark binding energies Ei are obtained by 

subtracting the internal energy of the basis state with the lowest energy, e.g. 

Ei = \- («i 3 + V24). 

A central element in the model is a phenomenological factor / appearing in the 
overlap of the basis states < A\B >= f/N c for SU(N C ). This factor is a function 
of the spatial coordinates of all four quarks, making the off-diagonal elements of 
V in eq. |2| four-body potentials, and attempts to take into account the decrease 
of overlap from the weak coupling limit, where < A\B >= 1/N C . Perturbation 
theory to 0(a 2 ) also produces the two-state model of eq. [I] with / = 1 Q. Several 
parameterizations for / have been suggested, a general form being 

/ = f ce - k AbsA~k P ^biP^ ^ 

Here 65 is the string tension, f c a normalization constant and k^, kp multiply 
the minimal area and perimeter bounded by the four quarks respectively. In this 
work either the normalization or the perimeter term is omitted from eq. £|[ 

This simple version of the "/-model" works for quarks in the corners of squares 
and (tilted) rectangles [(R) and (TR) in fig. |l|], but fails to predict some features 
of non-planar geometries, e.g. the doubly degenerate ground state energy of a 
regular tetrahedron [(T) with r = d]. In section [5] there will be introduced a 
generalization of the model that is capable of reproducing this degeneracy. 



3 Extrapolating to the Continuum 

3.1 Setting the scale 

Previously the scale in our simulations has been set by equating to the experi- 
mental continuum value the string tension bs = lim^oo F(r), where F(r) is the 
force between two static quarks, in lattice units. The string tension, however, was 
obtained by fitting the lattice parameterization of two-body potential to values 
from simulations at r/a = 2, . . . , 6 and not as r — > 00. 

Sommer has designed a popular new way to set the scale that uses only 
intermediate distances. First the force F(r/a) between two static quarks at 
separation r/a is calculated. By solving 

(r /a) 2 F(r /a) = c (5) 

with c = 1.65 for r /a, we get the equivalent of the continuum value r ~ 0.5 fm 
in lattice units. The constant on the right-hand side of eq. [5] has been chosen 
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to correspond to a distance scale where we have best information on the force 
between static quarks. According to the widely cited paper of Buchmiiller and 
Tye || various nonrelativistic effective potentials which successfully model heavy 
quarkonia agree in the radial region 0.1 to 1.0 fm and predict r.m.s. radii from 
0.2 to 1.5 fm for cc and bb systems. Even though these effective potentials are 
not the same as the QCD potential between static quarks, the distance scale 
where we have the best experimental evidence seems to be around 0.5 fm. A 
different estimate is given by Leeb et al. 0, who fit modified Cornell and Martin 
potentials to meson masses and claim that the known mesons determine the 
potential model-independently only between r = 0.7 fm and r = 1.8 fm. 

In the Cornell |J and Richardson || potential models the constant 1.65 in eq. 
|5| corresponds to r = 0.49 fm. On the other hand, the Martin model [HJ, with 



the strange quark counted as heavy, gives r = 0.44 fm and the modified Cornell 
and Martin potentials of ref. |7J result in 0.56 fm and 0.44 fm respectively (the 
published parameter values are incorrect [J). 

Choosing c = 2.44 as the scale setting constant makes the Richardson, modified 
Cornell and modified Martin potentials agree on r = 0.66 fm with the basic 
Cornell and Martin models giving r = 0.64 fm and 0.625 fm respectively. These 
models are clearly in better agreement at c = 2.44 than at c = 1.65. However, it 
is not excluded that the agreement may be accidental due to uncertainties in the 
models. 

To set the scale, a lattice parameterization of the two-quark potential was first 
fitted to values obtained from simulations at r/a = 2, . . . , 6 for each (3. The 
parameterization used was 

Vl{t) = - (^j + b s r + v , (6) 



where the on-axis lattice Coulomb potential is ITT 



1\ 3 r /a dhdfadfa e~ irk ' 



na 



rJ L J-Tr/a (2tt) 3 Ei=i sin (ha/2) 

„ 71 y cos(r gi ) 

L' s V£?=i sin2 (W2) 1 ) 

7T 7T(2L' - 1) 

m = — - — . (8) 

Here V is twice the number of spatial sites along one axis. The fit results are 
presented in table [1] for data from tilted rectangle and square geometry runs 
performed in Helsinki, and (3 = 2.5 square runs from Wuppertal. 



Personal communication with H. Leeb. 
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8 


g 


b?a 2 


Vq(1 


y 2 

A 


2.3 


299(14) 


1278(20) 


588(11) 


2.8 


2.35 


255(12) 


0985(15) 


557(9) 


2.6 


2 35fTR) 

— • '.J KJ \ X J- \j J 


D 260f5) 


0978(7) 


562(4) 


1 


2.4 


0.248(7) 


0.0709(10) 


0.554(6) 


0.1 


2.4(TR) 


0.238(9) 


0.0716(13) 


0.546(7) 


0.18 


2.45 


0.244(5) 


0.0494(7) 


0.551(4) 


0.04 


2.45(TR) 


0.238(5) 


0.0507(6) 


0.545(4) 


0.05 


2.5 


0.226(5) 


0.0373(7) 


0.531(4) 


0.55 


2.5(TR) 


0.233(4) 


0.0367(5) 


0.534(3) 


0.53 


2.5[W] 


0.234(4) 


0.0371(4) 


0.536(3) 


0.11 


2.55 


0.223(2) 


0.0271(3) 


0.522(2) 


0.03 


2.55(TR) 


0.224(5) 


0.0268(6) 


0.523(4) 


0.13 



Table 1: Parameters of the lattice two-quark potential. [W] refers to data from 
Wuppertal, (TR) to tilted rectangles. 



The force was calculated from this parameterization by the finite difference 



F W = V{r) - V / ~ d \ (9) 



using d = a. Here 



where 



ri -i[G{r-d)-G{r)Y (10) 

G{r) = - (-) (11) 

and Ti has been defined to remove 0[(a/r) 2 ] lattice artefacts from the argument 
of F, and to make F(ri) a tree-level improved observable ||. 

After getting the force at points rj/a, the expression (ri/a) 2 F(rj/a) was interpo- 
lated to get r /a corresponding to eq. §. Results for different /3's are presented in 
table ||, with lattice spacing a corresponding to r = 0.49 fm and a 11 to rj 1 = 0.66 
fm. For comparison, a value of a bs obtained from equating the string tension to 
the somewhat arbitrary continuum value \fb~s = 0.44 GeV is presented. In all 
cases the a 11 consistently agrees with a bs , whereas a disagrees, a bs /a being p& 1.09 
at each 3. Choosing ^fbs = 0.478(4) GeV would move a bs to the same value as 
a. A similar disagreement with the ratio of spacings being « 1.04 was found by 
the SESAM collaboration working in SU(3) and using c = 1.65 in eq. [ |T2|j . 
However, one should keep in mind that the quenched SU(2) string tension does 
not have to equal a phenomenological value. 
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As using c = 2.44 in the Sommer scheme produces better agreement with various 
continuum potential models and leads to agreement of a 11 and a bs at continuum 
string tension value b$ = 0.44 GeV, it seems that c = 1.65 underestimates the 
lattice spacing. Using c = 2.44, corresponding to tq « 0.66 fm, is in our case 
apparently a better choice than c = 1.65 and was chosen for this work, although 
the possibility that the difference is accidental cannot be excluded. 

The (3 = 2.5[W] values are calculated from potential data by a Wuppertal group 
doing simulations similar to ours fL3 |. The scales given by our and their data 
using the same analysis method agree well at (3 — 2.5 , but the scales obtained 
from the published b s values [0.0826(14) fm vs. 0.0864(5) fm] using the same 
data do not agree, which hints of the sensitivity of any intermediate distance 
determination of the string tension to the details of the fitting procedure. The 
errors in table |] were estimated by adding or subtracting the one-sigma errors of 
the parameters of the lattice two-body potential, so systematic errors from e.g. 
the choice of r are not included. 



p 


r /a 


a [fm] 


r o/ a 


a 11 [fm] 


a bs [fm] 


2.3 


3.25(4) 


0.1508(20) 


4.09(5) 


0.1613(18) 


0.1603(13) 


2.35 


3.76(2) 


0.1302(6) 


4.71(5) 


0.1401(14) 


0.1408(11) 


2.35(TR) 


3.77(2) 


0.1299(7) 


4.72(3) 


0.1397(6) 


0.1402(5) 


2.4 


4.45(5) 


0.1101(10) 


5.56(5) 


0.1186(10) 


0.1194(9) 


2.4(TR) 


4.44(2) 


0.1103(5) 


5.55(2) 


0.1190(4) 


0.1200(11) 


2.45 


5.34(4) 


0.0918(7) 


6.62(5) 


0.0997(8) 


0.0997(7) 


2.45(TR) 


5.27(5) 


0.0929(7) 


6.59(5) 


0.1002(7) 


0.1010(6) 


2.5 


6.18(5) 


0.0793(5) 


7.71(8) 


0.0856(9) 


0.0866(9) 


2.5(TR) 


6.21(5) 


0.0789(5) 


7.75(6) 


0.0851(6) 


0.0859(6) 


2.5[W] 


6.18(3) 


0.0793(3) 


7.71(5) 


0.0856(6) 


0.0864(5) 


2.55 


7.26(5) 


0.0675(5) 


9.04(6) 


0.0730(5) 


0.0741(7) 


2.55(TR) 


7.29(10) 


0.0673(9) 


9.08(12) 


0.0727(10) 


0.0734(8) 



Table 2: Values of r Q /a and a for each f3. r^/a and a 11 are calculated usin^ 
c = 2.44 in eq.|. 



The values of r obtained agree with the relation r = J(c — e)/bs using c from 
eq. |5] and e, bs from the fit of the two-body parameterization of eq. || 

A criterion for scaling of the two-body potential is equality of the parameter e 
for fits in the same physical distance range. In our case the fit range is shorter 
at higher /3's, and the values in table [I| support scaling at (3 = 2.5 and (3 = 2.55. 

In an attempt to take into account non-perturbative effects we can fit the values 
of a using the perturbative three-loop relation between a and (3 but with two 
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extra terms, 



a = -r e-»=*> (— p) 2 "o 1 + ^ ^^ 2N C /P + A 2 /(3* + A 3 /^ + 0(1/(3*] 

A-LAT P { ^Po 

(12) 

Here, as usual, (3q = y ^2 and (3\ = y(i^2-) 2 , while the recently calculated 



$ = ( I ^)3(_ 366 .2 + - ^3) g. Using Michael's result from ref. |5 

let us fix \ZI)s/Alat = 31.9. Fitting to all the values of a 11 in table || results 
in A 2 = -6.72(73), A 3 = 25.8(1.8) with X 2 = 0.9. The x 2 value is per d.o.f. as 



elsewhere in this paper. A modification of eq. [12] with no two to fourth order 
terms in 1/(3 but only a fifth order term with coefficient A 5 = 55.7(3) gives also 
a low x 2 = 

The /?-function b = d (3 /din a obtained using the former fit is presented in table 
|3| with comparisons to three-loop predictions and also estimates by Engels et al. 
using a thermodynamic approach Jl6| . The errors in our estimates were obtained 
by fitting to the values of a instead of a 11 in table ^ and using the form of eq. F2 
with only fifth order coefficient A 5 fitted in order to account for systematic errors 
from the choice of c in eq. [5] and the assumed functional dependence in eq. |T^. 

The agreement with Engels et al. is good except at the lowest (3 - the decrease in 
b with decreasing (3 at j3 < 2.37 is not reproduced by our approach. An estimate 
of b = —0.35(2) at (3 = 2.4 obtained using energy sum-rules [H] seems somewhat 
low when compared to the values in table |[ 

At (3 — 2.5 a fit of the two-loop analogue of eq. [L2] gives a value 78% of the 
perturbative two-loop prediction, which is modestly improved to 80% for the 
corresponding three-loop expressions. At these values of the bare coupling, addi- 
tional terms in the perturbative series are thus not likely to lead to a significantly 
better agreement with PT. 



p 


b from eq. |12| 


ref. 16] 


3-loop PT 


2.3 


-0.298(9) 


-0.3393 


-0.3898 


2.35 


-0.301(6) 


-0.3055 


-0.3896 


2.4 


-0.305(6) 


-0.3018 


-0.3893 


2.45 


-0.308(3) 


-0.3057 


-0.3891 


2.5 


-0.312(2) 


-0.3115 


-0.3889 


2.55 


-0.315(4) 


-0.3183 


-0.3887 



Table 3: Comparison between values of b = 8(3/ d In a from different methods. 
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3.2 Extrapolating two- body potentials 



To extrapolate to the continuum, we need values of energies at different /3's but 
corresponding to the same physical size. Two-body potentials were interpolated 
from the values given by the lattice parameterization. Figure |2] shows an extrap- 
olation of some of the two-body potentials involved in tilted rectangles. Others 



(from a total of 29) look similar - more pictures can be seen in ref. Jl8| . 

Linear (quadratic) fits with all data points included give x 2 values from 15 to 35 
(1.8 to 3.5), while cutting the (3 = 2.35 point off improves these to from 8 to 17 
(0.5 to 1.5). Continuum energies given by the latter quadratic extrapolations are 
6-8% higher than by the quadratic extrapolations including the (3 = 2.35 data. 

The 29 quadratically extrapolated continuum potentials can be parameterized 
using eq. ||with the lattice coulomb term replaced by e/r and e = 0.21(11), b$ = 
5.11(80) fm -2 and vq = 9.38(11) fm _1 (x 2 = 0.4). The value of the string tension 
corresponds to \/bs = 0.446(37) GeV, which agrees well with the continuum value 
used to obtain the a bs in table 0. 

The lattice artefact at each f3, i.e. difference from quadratically extrapolated 
continuum energy, has a linear dependence on the inverse separation of the quarks 
in lattice units squared (a/R) 2 . The slopes of this dependence at each (3 are within 
errors of each other, the average being 0.9(2). 



3.3 Extrapolating the binding energies 

The continuum extrapolations were performed for the square and tilted rectangle 
geometries, as the simple /-model works for them and allows the interpolation 
of binding energies to distances which are non-integer multiples of the lattice 
spacing. The ground state energy is given by 

E Q = {2v sl - 2v s2 + f 2 (v s2 - v d ) (13) 
+ yj[2v a2 - 2v sl - P(v s2 - v d )f - P(-A + P)(v d - v s2 y} 

-2 + 0.5/ 2 

where v s i, v s2 , v d denote the two-body potentials along the shorter and longer side 
and the diagonal of the rectangle, respectively. 

The two-parameter model of eq. § (with f c , or fc^, k P ) was fitted to the 
energies of squares and tilted rectangles separately. Result of extrapolating from 
these values is presented in figure § for tilted rectangles. Similar interpolations 
and extrapolations were performed for squares with sizes 2a, . . . , 6a at (3 — 2.5. 
Slopes for 11 of the 13 extrapolated energies were consistent with zero, supporting 
scaling of the binding energies at f3 > 2.35. Linear fits have low x 2 values and are 
more reliable than quadratic fits, which introduce an unnecessary extra parameter 
and lead to very large errors on the continuum values. 
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0.000 0.010 0.020 0.030 0.040 0.050 

(a/g 2 

Figure 2: r V vs. (a/r ) 2 for some two-body potentials, dimensions x, y, z shown 
in lattice units at f3 — 2.5 
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Figure 3: r E vs. (a/r ) 2 for tilted rectangles with dimensions in lattice units 
shown at f3 — 2.5 
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4 The /-model in the continuum 



4.1 Extrapolation of the parameters of / 

The parameters of the factor / were extrapolated from values given by a fit to 
the combined energies of squares and tilted rectangles, as shown in table |] for 
parameterizations with either the normalization or perimeter term omitted from 
eq. |. 



p 


fc 


k A 


x 2 


k A 


kp 


x 2 


2.55 


0.974(17) 


0.77(6) 


0.9 


0.71(10) 


0.020(16) 


1.0 


2.5 


0.992(17) 


0.76(4) 


0.3 


0.73(8) 


0.008(13) 


0.3 


2.45 


0.925(10) 


0.64(3) 


1.0 


0.49(4) 


0.057(8) 


1.2 


2.4 


0.864(15) 


0.57(2) 


1.25 


0.38(4) 


0.087(10) 


1.25 


2.35 


0.843(8) 


0.57(2) 


1.4 


0.36(3) 


0.101(5) 


1.1 



Table 4: Two-parameter /-models fitted to the square and tilted rectangle data. 



Typical values of / are quite constant at different /3's, as shown in fig. |], having 
slopes from —3.0(4) to 2.5(3). These slopes are a reflection of the behaviour of 
the two-body potentials used to calculate the /'s. The parameters f c , k A and kp 
are extrapolated in figure ||. 

The lattice artefacts of masses are known to behave roughly as a 2 , but no the- 
oretical justification exists for such behaviour of the artefacts of /. Therefore 
extrapolation results assuming this dependence should be given less weight than 
fits of the model to continuum extrapolated energies presented below. These two 
approaches agree for values of / in the continuum, while for the parameters of / 
the quadratic extrapolations roughly agree with values obtained from continuum 
fits of the model as can be seen in table [7]. 

4.2 Continuum fit 

The fit data for all linear extrapolations is presented in table |5| and for all 
quadratic extrapolations in table [(]. 

Fit and parameter extrapolation results can be seen in table [?|. Fits were also 
performed with quadratic extrapolations for two-body potentials and linear for 
four-bodies and are denoted "linear /quadratic" in the table. The parameter 
extrapolations have two \ 2 values, one for each of the extrapolated parameters. 
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Figure 4: A continuum extrapolation of /'s for squares with length of a side from 
2a (top) to 6a (bottom) at j3 — 2.5. 
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Figure 5: f c , &u and hp continuum extrapolations. 
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Dei [fill *1 


7; t ,o[firi * 


? } A f f rn ^ 1 


2,2,0 


0.92(6) 


-0 515(25) 


7.44(34) 




8 26(30) 


3 2 2 


0.88(6) 


-0 379(15) 


8 34(39) 


8 432(40) 


9.14(40) 


3,3,0 


0.77(2) 


-0 439? 101 


8 38(39) 




9 17(39) 


3 3 1 


84(4) 


-0 358 (1 1 ) 


8 43(40") 


8 55(40") 


9 28(40) 


4 3 2 


0.69(3) 


-0.261(12) 


8.84(40) 


9.07(40) 


9.84(40) 


4,4,0 


0.62(2) 


-0.430(9) 


8.98(31) 




9.89(31) 


4,4,1 


0.65(2) 


-0.376(15) 


9.08(40) 


9.15(40) 


10.02(40) 


5,4,3 


0.48(2) 


-0.406(19) 


9.63(40) 


9.63(40) 


10.68(40) 


5,5,0 


0.46(2) 


-0.388(14) 


9.55(30) 




10.59(31) 


5,5,1 


0.48(2) 


-0.354(20) 


9.63(40) 


9.69(40) 


10.71(40) 


6,5,3 


0.35(2) 


-0.281(27) 


10.07(40) 


10.15(40) 


11.30(40) 


6,6,0 


0.33(2) 


-0.334(19) 


10.06(31) 




11.25(32) 


6,6,1 


0.32(2) 


-0.296(25) 


10.15(40) 


10.19(40) 


11.38(40) 



Table 5: Continuum fit data with linear extrapolations 



All fits and extrapolations predict f c to be one or very slightly above. If it is 
set to one the fit to linear data and linear extrapolation give a fc^ somewhat 
below one, while the linear/quadratic fit has k^ approximately one and the fit 
to quadratic data has a value slightly above unity. Since the linear/quadratic 
fits have low x 2 values and are most reliable for two-body potentials and four- 
quark binding energies, our best estimates for the continuum values are f c = 1, 
kji = 1.0(1) and kp = 0. The perimeter term should become unimportant in the 



continuum limit [D5|, which is indeed observed. Thus the only parameter with 
physical significance is k^- 

Continuum fits were also performed for squares extrapolated using the original 
value c = 1.65 in Sommer's equation [5]. The resulting values of f c and k^ for 
all linear and quadratic fits were within errors with fits to extrapolations using 
c = 2.44. As another test linear /quadratic and quadratic fits were performed with 
two-body potentials extrapolated with the (3 = 2.35 point cut off, with similar 
results but higher \ 2 values. 

In the strong coupling limit a factor exp (—bsA) appears in the diagonal elements 
of the Wilson loop matrix, A being the minimal surface bounded by straight lines 
connecting the quarks. The factor is closely related to / and originates from the 
sum of spatial plaquettes tiling the transition surface between two basis states. 
When moving to weaker couplings there exist correlated flux configurations for 
which the transition area is much smaller than A, which has been expected to lead 



to larger mixing among basis states [20]. If this is the case, it is not reflected in 



our best estimate of a kA about one. A larger transition area could be explained 
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x / o, 11 n, z 1 n, 


J K^OJ 


_L/(J |±1±1 j 


7) n1 [frn M 


7? offm ^ 

u s2 [mi 




2,2,0 


90(16) 


-0.40(11) 


9 29(17) 




9 95(18) 


3 2 2 


86(8) 


-0 46(13) 


9 70(18) 


9 73(18) 


10 52(20) 


3,3,0 


0.75(11) 


-0 43(10) 


10 09(19) 




10 88(21) 


3 3 1 


68(5) 


-0 29(14") 

VJ . Zj I It: f 


9 73(18) 


9 89(19) 


10 65(20) 


4 3 2 


0.55(6) 


-0.20(11) 


10.22(19) 


10.45(20) 


11.22(22) 


4,4,0 


0.56(9) 


-0.40(10) 


10.68(20) 




11.59(23) 


4,4,1 


0.52(5) 


-0.30(15) 


10.45(20) 


10.53(20) 


11.40(23) 


5,4,3 


0.37(4) 


-0.33(19) 


11.00(22) 


11.00(22) 


12.07(25) 


5,5,0 


0.35(7) 


-0.32(11) 


11.24(22) 




12.31(26) 


5,5,1 


0.36(4) 


-0.27(18) 


11.00(22) 


11.07(22) 


12.10(25) 


6,5,3 


0.27(4) 


-0.23(18) 


11.46(23) 


11.52(23) 


12.69(28) 


6,6,0 


0.19(5) 


-0.22(11) 


11.76(24) 




12.99(29) 


6,6,1 


0.21(4) 


-0.20(29) 


11.52(23) 


11.56(23) 


12.77(28) 



Table 6: Continuum fit data with quadratic extrapolations 



by the finite width of actual flux tubes instead of simply lines as in strong coupling 
approximation. 

5 How can the model be developed? 

The failure of the simple /-model of eqs. p] and |^ to predict the energies of 
the tetrahedral geometry has been proposed to be due to a dependence of the 
four-quark energy on the first excited state(s) of the gluonic field in the static 
two-quark potential []]]. A higher lying basis state would make the ground state 
attractive with increasing attraction when the size of the system gets larger, since 
the excited state gets closer to the ground state energy. Such a trend of increasing 
attraction is observed in the tetrahedral energies unlike in the case of squares, 
whose energies get smaller with increasing size. 

The lowest-lying gluonic excitation has non-zero angular momentum about the 
interquark axis, transforming in the E u representation of the cubic symmetry 
group Dihi while the second lowest excitation has the same A\ g symmetry as the 
ground state. In our calculation the two-body paths used as a variational basis 
differ in the amount of fuzzing (a.k.a. smearing, blocking) applied to the lattice. 
Fuzzing only changes the overlap with the excited states (to be denoted with a 
prime) possessing the A\ g symmetry and does not bring in any of the lowest E u 
excited state. The question remains which one of these excited states has more 
effect on the tetrahedral energy. 
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f, 

J c 




k P 


y 2 

A 


all linp^r 

(111 ±±J.J.V_vC1j-1- 


1.08(3) 


1.04(5) 


o 


0.8 


all linear, f c = l 


1 


0.94(3) 


o 


1.3 


linear extrapol. 


1.04(3) 


0.83(6) 





2.6/3.3 


linear extrapol. (II) 


1 


0.85(9) 


-0.02(2) 


2.3/2.6 


linear/ quadratic 


1.04(4) 


1.04(6) 





0.7 


linear/quadratic, f c — 1 


1 


0.98(3) 





0.8 


all quadratic 


1.08(8) 


1.39(13) 





0.4 


all quadratic, f c — 1 


1 


1.27(6) 





0.5 


quadratic extrapol. 


1.10(5) 


1.09(11) 





2.7/0.4 


quadratic extrapol. (II) 


1 


1.2(2) 


-0.08(4) 


0.8/1.6 



Table 7: Continuum fit and parameter extrapolation results. 



The dependence of the four-quark energies on the A' lg excited states can be 
investigated using the fact that an increasing fuzzing level reduces the overlaps 
Cj, i > 1 of these excited states in the energy eigenstate expansion of the Wilson 
loop 

< W(R, T) >= C n {R)e- y ^ R)T . (14) 

n 

The overlaps C n (R) > have a normalization condition 

£c„(i?) = i. 

n 

Increasing the fuzzing level used in the calculation of four-body energies should 
worsen the convergence of any four-body energy with a significant dependence on 
the excited state(s) with the symmetry of the ground state. Our group has been 
using fuzzing level 20 with the factor multiplying the link to be fuzzed c = 4, 
while a group in Wuppertal doing related work has chosen 150 with c = 2 to 
maximize the ground state overlap. At (3 — 2.5 they obtain Cq values from 0.81 
to 0.98 on a 16 4 lattice. 

We made runs at — 2.5 on a 24 3 x 32 lattice with three different fuzzing levels 
used to calculate the four-quark energy of a regular tetrahedron. Ground state 
overlaps Co of two-quark paths at different separations r/a,d/a and at different 
fuzzing levels are shown in table |[ The four-body operators are constructed from 
the same paths, and the relative quality of convergence of the binding energy E± 
as a series in T is also shown. 

The table shows that best convergences of the four-quark energies are obtained 
when the two-body paths have practically no overlap with gluonic excitations. 
Therefore the only way excitations can contribute is through overlaps of two- 
quark paths with excited states of gluon fields between other quark pairs. This 
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contribution is likely to be significant; e.g. the Isgur-Paton string model pi 



with N=2 predicts overlaps between Is, lp (corresponding to A\ g , E u ) and Is, 2s 
(corresponding to A\ g , A' Xg ) states to be ~ 0.4 when the centres of two parallel 
fluxtubes are separated by ~ 0.4 fm and ~ 0.7 fm respectively. This can be 
seen from fig. |6]. The rapid worsening of convergence in table |8| with decreasing 
ground-state overlap and increasing overlap with A' lg excitations suggests that the 
lowest-lying gluonic excitation with E u symmetry has a more important effect on 
the four-quark energies than the higher lying A' lg excitations. 





FL=150 


FL=20 


FL=0 


r/a, d/a 


C % 


C % 


C % 


2,2 


N/A 


100.0(1) 


98.2(1) 




interm. 


best 


worst 


3,3 


98.7(1) 


100.0(1) 


96.7(1) 




interm. 


best 


worst 


4,4 


98.2(2) 


99.7(1) 


94.3(1) 




"best" 


"best" 


worst 


5,5 


97.9(2) 


99.5(1) 


91.1(1) 




interm. 


best 


worst 


6,6 


98.0(2) 


99.1(1) 


85.3(1) 




interm. 


best 


worst 


7,7 


99.1(4) 


98.0(1) 


N/A 




interm. 


best 


worst 



Table 8: Ground state overlaps and E4 convergence at different fuzzing levels 
(FL) and interquark distances. 



The simple /-model of eqs. [I] and |2| can now be extended by using (in addition to 
the A, B, C shown in figure P three additional basis states describing the same 
quark partitions but involving the first excited state of the two-body potential. 
The resulting 6x6 matrices contain two new gluon field overlap factors f a and 
f c similar to /, which however measure the overlap between an excited and 
ground state basis state and two excited basis states respectively. Solving the 
determinant analogous to eq. [I] gives for a regular tetrahedron two degenerate and 
two non-degenerate energy eigenvalues. The degenerate negative ground state 
energy of a regular tetrahedron would correspond to the lower of the degenerate 
eigenvalues, while it is harder for the model to predict the third eigenenergy from 
simulations since it is dominated by gluonic excitations as argued in ref. ]l[. 

To estimate these eigenvalues, a parameterization for / from table |] can be used 
with the minimal transition areas for regular tetrahedra calculated in ref. |p2 |. 
A parameterization for the energy of the first excited two-body state is found 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
fluxtube separation in fm 



Figure 6: Overlaps of fluxtube wavefunctions as functions of fluxtube separation 
in the Isgur-Paton model with N=2. 



from ref. ||17|| . Because all the quarks in a regular tetrahedron have the same 
distance from each other, all the excited basis states have the same energy. Using 
these values, the simulation results for regular tetrahedra and a self consistency 
argument allows estimation of other parameters of the model. 

The values obtained give the lower degenerate eigenvalue as the lowest eigenvalue 
when f c > 1 because of a non-degenerate eigenvalue being unstable around f c = 
1. This unstable eigenvalue is probably stabilized by a minor contribution from 
the A' lg excited state, requiring treatment with 9 basis states. 

The values for the new overlap factor f a have a similar linear behaviour with the 
size of the system as is observed in figure || in the overlap of ls-lp Isgur-Paton 
fluxtubes with separation in the same distance range, supporting reliability of 
the estimation scheme. 

Work is in progress on applying this extended model to geometries more compli- 
cated than the regular tetrahedron . 
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6 Conclusions 



The conclusions in this work can be summarized as follows: 

1. Choosing c = 2.44 (corresponding to r ~ 0.66 fm) in eq. |5| leads to better 
agreement of various continuum potential models with each other, and in 
our case gives lattice spacings that agree with those obtained from the string 
tension. This value of c thus seems a better choice than the original choice 
c = 1.65 (t~o ~ 0.49 fm) by Sommer, although it cannot be excluded that 
the difference is accidental due to uncertainties in the potential models and 
the quenched SU(2) string tension. 

2. The lattice spacings a at the six values of (3 used in simulations give the f3- 
function b = in agreement with ref. |16j when /3 > 2.35. At f3 — 2.5 the 



measured value accounts for 80% of the three-loop perturbative prediction, 
presenting a very modest improvement of 2% from the two-loop case. 

3. A quadratic extrapolation is preferred for two-body potentials, whose con- 
tinuum parameterization is given in section [3.2| . Four-quark binding ener- 
gies are scaling from (3 > 2.35 and their values do not change significantly 
in the range of (3 values simulated. In practice the simulation results for 
the binding energies can thus be used directly, whereas it is recommended 
to use the continuum parameterization for two-body potentials to avoid 
introducing lattice artefacts into a four-quark energy model. 

4. Parameter extrapolations and continuum fits of the simple /-model show 
that in the continuum, the normalization of the gluon field overlap factor / 
can be safely set to one and the perimeter term to zero, leaving the constant 
kj{ multiplying the area as the only physically relevant parameter. Its value 



is about one and not C 1 as predicted for the transition area in ref. [ 20 1 



5. In our simulations, the effect on the four-quark energies of excited states 
of the gluon field between two quarks comes from the overlap of two-body 
paths with excited two-body fields of other quark pairings. An extended 
/-model using the excited states of the two-body potentials can reproduce 
the negative, degenerate ground state energy of a regular tetrahedron unlike 
the simple model based on ground state potentials. 

Other possible extensions of the simple /-model include instanton effects (ref. 
24]]) and four-body interactions in the strong coupling limit (ref. |2f 



A study using microscopical flux distributions would be helpful in determining 
the relation of the shape of two- and four-quark fluxtubes and their overlaps to 
the parameterization of /. Our next project attempts this by using sum-rules 
similar to those derived by C. Michael EH . 
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